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Abstract. We study Ising ferrimagnets on square lattices with antiferromagnetic 
O , exchange couplings between spins of values S=l/2 and S=l on neighbouring sites, 

couplings between S=l spins at next-nearest neighbour sites of the lattice, and 
a single-site anisotropy term for the S=l spins. Using mainly ground state 
considerations and extensive Monte Carlo simulations, we investigate various aspects of 
the phase diagram, including compensation points, critical properties, and temperature 
dependent anomalies. In contrast to previous belief, the next-nearest neighbour 
couplings, when being of antiferromagnetic type, may lead to compensation points. 
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1. Introduction 



Mixed spin ferrimagnetic Ising models have been studied for some time, with a renewed 
interest, especially, in connection with 'compensation points', see, for instance, fi- 
ll]. These points occur at temperatures, below the critical temperature, at which the 
sublattice magnetizations cancel exactly, giving zero total moment. As the temperature 
is tuned through such a point the total magnetization changes sign, which may be used 
in technological applications, most notably in magnetic recording. 

One of the simplest such models consists of classical Ising spins S=l/2 and S=l 
on a square or simple cubic lattice with the spins of different type being located on 
neighbouring sites of the lattice. Its Hamiltonian may be written in the form 

H = Jx 5>,S, + I> £ £/» (1) 

(hj) iefi 

where J\ > denotes the antiferromagnetic coupling between spins Oj = ±1/2 on 
the sites of sublattice 'A', and neighbouring spins Sj = 1, 0, —1 on sites forming the 
sublattice 'B'. Then, spins on each sublattice tend to order ferromagnetically, with 
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opposite sign for the two types of spins. D is the strength of a single-site anisotropy 
(or crystal-field) term acting only on the S=l spins of sublattice B. 

One may also choose <7j = ±1 rather than ±1/2. The change has to be taken 
into account when calculating sublattice magnetizations and defining the compensation 
point. Otherwise, the modified convention simply amounts to a rescaling of the exchange 
coupling [51 [10] . 

The mixed spin Ising model, eq. (1), is known, see, e.g., [U El [10], to lead to 
compensation points for simple cubic lattices, in accordance with mean-field theory [2] , 
but not for square lattices, in contrast to mean-field theory. As has been observed 
already some years ago by Buendia and Novotny j5j [12], compensation points may 
occur on square lattices, when adding a (ferromagnetic) coupling between next-nearest 
neighbouring (nnn) spins on the A sublattice. On the other hand, the authors did not 
find any evidence for compensation points, when considering nnn couplings for B spins, 
instead of the ones for A spins. 

In the following article we shall challenge the latter suggestion, which seems to have 
been taken for granted by others, see, e.g., [7J. Indeed, based on extensive Monte Carlo 
(MC) simulations, we shall present clear evidence for compensation points due to nnn 
antiferromagnetic interactions between B, or S=l, spins on the square lattice. In fact, 
the effect is not easy to identify. 

The outline of the article is as follows. In section 2, to set the scene for the following 
parts, we define the model with nnn couplings between S=l spins, discuss its ground 
states, and describe details of the simulations. The following section deals with the 
compensation points. In section 4, MC results on critical properties and temperature 
dependent anomalies of the model will be presented. A brief summary is given in the 
final section. 

2. Model, ground states, and simulations 

We shall study the following mixed spin Ising model on a square lattice 

n = j 1 J2^ i s j + Dj2s j 2 -J2 E ( 2 ) 

with antiferromagnetic nearest neighbour couplings, J\ > between A spins, cr, = ±1, 
and B spins, Sj = 0, ±1, a crystal-field term of strength D acting on the B spins, and 
nnn couplings, J2, between B spins. Note that we set A spins equal to ±1, in accordance 
with previous work J5J[T0]. To identify possible compensation points, care is then needed 
in defining the magnetization of the sublattice A, see above and below. 

To determine the ground states of the model in the (D/ J l5 J%j J\) plane, we first 
select, like before [HI [121 EB], the structures which are stable among the ones described 
by 2x2 cells of spins on the square lattice. One readily observes the possibility of 
degenerate ground states, eventually with indefinitely large unit cells. Finally, we check 
the tentative ground state structures by monitoring spin configurations at very low 
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temperatures obtained from careful cooling MC runs for lattices of various sizes. The 
resulting ground state phase diagram is depicted in figure 1. 

The phase diagram at vanishing temperature, T = 0, comprises four structures. The 
antiferromagnetic (AF) structure is stable for ferromagnetic or weakly antiferromagnetic 
nnn couplings, J2, and negative or sufficiently small positive values of the crystal field, 
D. The spins on each of the two sublattices order ferromagnetically, ±1, with opposite 
sign on the sublattices A and B. 

In two of the remaining other three ground states, S=l spins may be in the state 0, 
due to the single-site anisotropy term, D. Actually, for sufficiently large values of D, all 
spins of the sublattice B are in state 0. We abbreviate the resulting structure as 'OIF. 
The structure is highly degenerate, with each A spin being either —1 or 1, leading to a 
2 Ar - A -fold degeneracy for Na sites on the sublattice A. In contrast, in the 'OF structure, 
see figure 1, only half of the S=l spins are in the state 0. They form, on the sublattice 
B, in the standard Wood's notation for overlayers [Hj, a periodic c(2 x 2) superlattice. 
The other half of the spins on the B sites are ferromagnetically ordered, ±1. The spins 
on the sublattice A are ferromagnetically ordered as well, with opposite sign, compared 
to that of the B spins. 




D/Jj 

Figure 1. Ground state phase diagram of the mixed spin model, eq. (2), on a square 
lattice, with solid lines separating the different phases. The dashed line in the 01 phase 
refers to equation (7), see text. 

In the fourth ground state structure, all spins on the sublattice B are 
antiferromagnetically aligned, due to the dominant antiferromagnetic nnn coupling, J 2 . 
The structure may be abbreviated as B-AF. Each spin on the sublattice A may be either 
— 1 or 1, leading, as for Oil, to a large degeneracy. 

Note that there may be additional degeneracies at borderlines between different 
ground state structures [ID] . 

To determine thermal properties of the model, (2), we mainly perform MC 
simulations, using the Metropolis algorithm with single-spin flips, providing, indeed, the 
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required accuracy, so that there is no need to apply other techniques like cluster-updates 
or the Wang-Landau approach [15] . We consider lattices with LxL sites, employing full 
periodic boundary conditions. L ranges from 4 to 120, to study finite-size effects. 
Typically, we do runs of 5 x 10 6 to 10 7 Monte Carlo steps per spin, where averages and 
error bars may be obtained from evaluating a few of such runs, using different random 
numbers. These quite long runs lead to reliable data, as before |10| . The estimated 
errors are usually smaller than the sizes of the symbols in the figures, and they are 
shown there only in a few cases. 

We record the energy per site, E, the specific heat, C, both from the energy 
fluctuations and from differentiating E with respect to the temperature, and the absolute 
values of the sublattice magnetizations of the two sublattices 

I^a|=<|E^I>/( 2 ( L V2)) (3) 

A 

and 

\m B \=<\^S J \>/(L 2 /2) (4) 

B 

as well as the absolute value of the staggered magnetization of sublattice B, describing 
the ordering for the B-AF structure, 

KI=<IE^-E^I>/(^ 2 /2) (5) 

B+ B~ 

with B + '~ denoting an obvious bipartition of sublattice B. The brackets <> denote 
the thermal average. Note the factor of 1/2 in the definition of taking into 

account the correct length of the S=l/2 spins, so that |m^(T = 0)| = 1/2 for the 
ferromagnetic ground state of sublattice A. In addition, the corresponding sublattice 
(staggered) susceptibilities, xa, Xb, and Xbi have been computed from the fluctuations 
of the (staggered) sublattice magnetizations. We also analyse the fourth-order cumulant 
of various (sublattice) order parameters, the Binder cumulant [IS] , defined by 

U = 1- < m 4 > /(3 < m 2 > 2 ) (6) 

with < m 2 > and < m 4 > being the second and fourth moment of (staggered) 
magnetizations of sublattice A or B. Finally, we monitor typical equilibrium Monte Carlo 
configurations, illustrating the microscopic behaviour of the system and providing, e.g., 
evidence for the ground states structures, as mentioned above. 

To test the accuracy of the simulations, we compared our MC data to those of 
previous accurate numerical work and to exact results by enumerating all possible 
configurations for small lattices with L = 4 (5j [10] ■ 
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3. Compensation points 

Compensation points occur at temperatures T comp , where both sublattice magnetizations 
cancel each other, with a vanishing total magnetization. In finite systems, as one is 
studying in MC simulations, a convenient and efficient way [5] to locate such points is to 
use the crossing condition |m^|(T comp ) = \mB\(T comp ). At low temperatures sufficiently 
far below the phase transition, finite-size effects are expected to be very weak, because 
there the compensation effect is not related to critical phenomena. 




Figure 2. Sublattice magnetizations \itia\ (open symbols, dotted lines) and \itib\ (full 
symbols, solid lines) at fixed D/J — 3.0 and varying J 2 / ' J\ = 0.25 (circles), 0.0 (triangles 
down), — O.f (squares), —0.24 (diamonds), and —0.26 (triangles left). Lattices of size 
60 2 are simulated. 

In the AF ground state, one has \7tib\{T = 0)=1, while |m^|(T = 0)= 1/2. Now, 
at non-zero temperatures, |m#|(T) may fall off quite rapidly due to antiferromagnetic 
couplings J2 and due to a relatively large single-site anisotropy term, D, which favours 
flips of S=l spins to state 0. Furthermore, at zero temperature, |j71b|(T = 0) drops from 
1 to 1/2, when passing through the borderline between the AF and 01 ground states. 
Accordingly, one might speculate that compensation points, at low temperatures in the 
AF phase, may show up close to the AF-0I borderline at T=0. In fact, our preliminary 
mean-field calculations seem to suggest that lines of compensation points may spring 
from this AF-0I borderline. 

However, in our MC simulations, investigating carefully several cases, we find no 
evidence for such compensation points in the AF phase. An example is depicted in figure 
2. There, D/J\ is fixed at 3.0, and Jij J\ is varied, from +0.25 to —0.26, so that the 
AF-0I border at J2/J1 = —0.25 is approached and crossed. One observes, that \w,b\ is 
always larger than |m^|, albeit the difference between the two sublattice magnetizations 
may get quite small when decreasing J 2 . MC data for lattices of fixed size, L = 60, are 
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displayed. The observation on the absence of compensation points in the AF phase also 
holds for other lattice sizes as well as for other values of J 2 j ' J\ and D/Ji. 




Figure 3. Sublattice magnetizations (open symbols, dotted lines) and ms| 

(full symbols, solid lines) at fixed J2/J1 = —0.2, and varying D/Ji= 3.0 (circles), 3.4 
(squares), 3.57 (diamonds), and 3.7 (triangles up), simulating lattices of size L= 60. 
At D/J\ =3.57, a compensation point shows up, T comp . 

On the other hand, compensation points will be argued below to occur in the 01 
phase, for J2/J1 > —1.0, arising at vanishing temperature from the line 

J 2 /J, = -2 + D/(2Ji) (7) 

which is depicted as the dashed line in the ground state phase diagram, figure 1. 

Let us first present numerical support for this claim, followed then by low 
temperature energy considerations backing it up. A numerical example is shown in 
figure 3, where, for lattices with 60 2 sites, the temperature dependences of the sublattice 
magnetizations, \m^\ and |?71b|, are shown at fixed nnn coupling, J 2 /J± = —0.2, and 
varying the strength of the single-site anisotropy term, D/Ji, from 3.0 to 3.7. Note 
that for this ratio of couplings, at T = 0, the AF-0I border is at 3.2, and the origin, 
at vanishing temperature, of the line of compensation points would be, according to 
equation (7), at D/J= 3.6. As discussed in the context of figure 2, one observes 
that \m,B\(T) is always larger than 1 771,4 ICO when being in the AF phase. However, 
in the 01 phase, here, at D/Ji= 3.57, there seems to be a compensation point, with 
kBTcomp/ J\ ~ 0.26, clearly below the phase transition: At lower temperatures, |ttib| is 
still larger than |m^|, with reverse ordering at T > T comp . By further increase of D, 
D / Ji= 3.7, I m^i I (T) is larger than |m B |(T), at all temperatures up to T c . 

To establish numerically the presence of compensation points, careful finite-size 
analyses are required, as illustrated in figure 4. Here, MC data for various lattices sizes, 
with L ranging from 20 to 120, are displayed, at J 2 j J\ = —0.2 and D / J\= 3.59, i.e. very 
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Figure 4. Sublattice magnetizations \itla\ (open symbols, dotted lines) and |ms| 
(full symbols, solid lines) at J2/J1 = —0.2 and D/J\= 3.59, with lattice sizes L= 20 
(squares), 40 (diamonds), 60 (triangles left), 80 (triangles down), and 120 (circles). 
The compensation point is located at ksT '/ ' J\ ss 0.16. 



close to the point, from which, according to (7), the line of compensation points may 
arise. In fact, a compensation point may be located at /csT comp /Ji w 0.16. Below that 
temperature, \m B \{T) supercedes \m,A\{T), with finite-size dependences being extremely 
small. Then, at T comp , a crossing of the two sublattice magnetizations occurs. The finite 
size effects are still very weak up to ksT/Ji ~ 0.25, providing clear evidence on the 
existence of the compensation point in the thermodynamic limit. 

Indeed, similar finite-size analyses allow one to locate the rather steeply rising line of 
compensation points, at fixed nnn couplings J 2 / J\ = —0.2. Results on the compensation 
points, T comp , are summarized in figure 5, together with estimates for the transition line, 
T c , to the paramagnetic phase. 

The critical line, T c , has been determined by monitoring the size dependent 
positions of maxima in the specific heat C, in the susceptibility of the sublattice A, 
and in the intersection points of the Binder cumulant, especially, of sublattice A. The 
intersections of the Binder cumulant U(L, T) for successive lattice sizes [16] turn out 
to have the smallest finite-size effects, thereby being most efficient in estimating the 
critical temperature. 

The line of compensation points seems to start at zero temperature at D/J\— 3.6 
for J%j J\ = —0.2, in agreement with equation (7). Furthermore, it extends only over 
quite a small region, with the compensation point coinciding with the critical point at 
about D/J l = 3.52 ± 0.02. 

Let us now turn to the low temperature considerations leading to (7). Specifically, 
we consider, for the 01 structure, the one-spin flip excitations on sublattice B. For 
concreteness, we analyse ground states with spins being in state —1 or on sublattice 
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Figure 5. Phase boundary to the paramagnetic phase, T c (solid line), and 
compensation points, T comp (dashed line), at J2/J1 = —0.2, varying Dj J\. 




Figure 6. Phase boundary to the paramagnetic phase, T c (solid line), and 
compensation points, T comp (dashed line), at J2/J1 = —0.8, varying D/J\. 



B, and in state 1 on sublattice A. Then there are four possible flips of S=l spins, 
namely flipping a spin from its state to the state (i) —1 or (ii) 1, or flipping a spin 
from its state — f to state (iii) or (iv) +1. Obviously, only flip (i) will increase the 
sublattice magnetization |ms| above its value at zero temperature, |me|(T = 0) = 1/2, 
while otherwise \m B \(T) will have a negative slope at low temperatures. Now, one 
may readily calculate the energies of all four elementary flips. One obtains that flip 
(i) costs the lowest energy either followed by flip (iii), when D/Ji < 2J 2 /J\ + 4, 
or followed by flip (ii), when J2/J1 > —1.0. Equation (7) is then obtained under 
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the assumption, that an initial increase of \7tib{T)\ is needed to have a compensation 
point, with |m_e| > |m^| at sufficiently low temperatures. Otherwise, in the 01 phase, 
\tua\ > \itib\ at all temperatures < T < T c , excluding a compensation point. 

We further checked numerically the hypothesis, (7), in determining lines of 
compensation points by fixing the nnn coupling J 2 J J\ not only at —0.2, but also at 
—0.4 and —0.8, varying the strength of the single-site anisotropy term, DjJ\. Indeed, 
we observe compensation points in the 01 structure, which seem to arise, at zero 
temperature, from D / J\= 3.2 at J2/J1 = —0.4 and from D / J\= 2.4 at J%j 1 J\ = —0.8, 
in agreement with (7). The latter case is displayed in figure 6. 

Note that the range of values of D/J±, in which there is a line of compensation 
points, shrinks appreciably as one decreases Jij J\- This may be seen by comparing 
figures 5 and 6. Thence, it becomes increasingly difficult to locate compensation points 
for lower values of J%j J\. 

4. Phase transitions and anomalies 

While mixed spin Ising models are of much interest because of the potential presence of 
compensations points, they may exhibit other intriguing thermal properties as well. 

We shall focus on two aspects, the characterization of the phase transitions 
associated with the various ground states, and anomalies in the temperature dependence, 
especially, of the sublattice magnetization and the specifc heat C. 

One expects that there is no phase transition associated with the highly degenerate 
Oil ground state for the following reasons: At zero temperature, all spins of sublattice 
B are in state 0, while each spin of sublattice A may be either —1 or 1. The single-site 
anisotropy term, favouring the spin state on sublattice B, does not support long-range 
order, in close analogy to the situation in the Blume-Capel model [17]. The spins on 
sublattice A act merely like a fluctuating random field, and may not lead to a phase 
transition neither. In fact, our MC data, for various quantities, like the specific heat C 
or the probability to encounter a spin in state 0, give no indication for singular thermal 
behaviour. 

In case of the other three ground states, AF, 01, and B-AF, we determine the 
universality class [TH [19] by analysing specific heat and susceptibilities. In particular, 
we monitor the size, L, dependence of that maximum of the specific heat, C max (L), which 
goes over into a singularity in the thermodynamic limit, L — > 00 (note that there may be 
other maxima as will be discussed below in the context of anomalies). Furthermore, we 
study the size dependence of the maxima in (sublattice) susceptibilities. For the AF and 
01 structures, we focus on XA,max(L). For the B-AF structure, we record the staggered 
susceptibility of the antiferromagnetically ordered sublattice B, Xb ma x(,L) ■ In all three 
cases, the critical behaviour is found to be consistent with having phase transitions in 
the universality class of the standard two-dimensional Ising model. 

Typical MC results on the size, L, dependence of maxima in the specific heat, 
C max , are depicted in figure 7. For all three types of ground states, AF, 01 as well 
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Figure 7. Height of the critical maximum of the specific heat, C max versus logarithm 
of the lattice size, InL, for L= 10, 20, 40, 60, and 80 at (a) (circles) J2/J1 = —0.3, 
D/J\ = 1.5, (b) (squares) J2/J1 = —1.0, D/J\— 1.5, and (c) (diamonds) J2/J1 = —1.5, 
D/Ji— 2.5. The three cases refer to the (a) AF, (b) 01, and (c) B-AF structures. 



as B-AF, one approaches to a good degree, already for lattices of moderate sizes, the 
form C max oc InL, being characteristic for the two-dimensional Ising universality class. 
Note that the prefactor in front of the logarithmic term depends strongly on J2/J1 and 
D/Ji. Certainly, such Ising-like critical behaviour may be expected in the AF case. In 
case of the B-AF phase, the spins on sublattice B order antiferromagnetically, with the 
A spins providing effectively a randomly fluctuating field, which may be argued to be 
irrelevant for the universality class. Finally, in the 01 case, the Ising-like criticality may 
be understood by the ferromagnetic order of the spins on sublattice A. 

The analysis of the susceptibilies confirms the findings on the specific heat. Typical 
examples, for the same model parameters as in figure 7, are shown in figure 8. The size 
dependent maxima in the susceptibilities are observed to follow the form Xmax oc L 7 ^ 4 , 
with rather small corrections, for all three cases, AF, 01, and B-AF. This behaviour 
may be quantified by calculating the slope between successive points in the doubly 
logarithmic plot of the MC data. The resulting effective local critical exponent of the 
susceptibility is near 7/4, even for moderate lattice sizes. Thence, in all three cases, 
criticality seems to belong to the two-dimensional Ising universality class. 

Care is needed when attempting to determine the universality class from the value 
of the critical Binder cumulant U*= U(T C ,L — > oo), because that value is known to 
depend on boundary conditions, shape, and anisotropy of the correlations (EH [201 EI]- 
Note that in the cases we studied, its value seems to be close to that of the standard two- 
dimensional Ising model with periodic boundary conditions for lattices of square shape, 
U* ~ 0.6107 [22]. However, the dependence, especially on anisotropy, may be very weak 
|20[ [21] . and a reliable analysis may require an exact or, at least, extremely accurate 
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Figure 8. Log-log plot for maxima of (staggered) sublattice susceptibilities Xmax 
versus lattice size L for the AF, 01, and B-AF cases, with the same parameters and 
notation as in figure 7. In case of AF and 01 structures, XA,max, in case of the B-AF 
structure, Xb max ^ s recorded. For comparison, the solid line shows Xmax °c L 7 / 4 . 



determination of the critical point. We therefore refrained from such an analysis here. 




Figure 9. Specific heat C as a function of temperature, fc^T/Ji, at J2/J1 = —0.2 and 
D/J = 3.4, simulating lattices of sizes L= 20 (circles), 40 (squares), 60 (diamonds), 
and 80 (triangles left). 

Let us now turn to discussing anomalies, exhibiting intriguing non-monotonic 
temperature dependences, which are not related to phase transitions. Especially, we 
monitor the magnetizations of sublattice B, |ms|, the specific heat C, and the thermally 
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averaged occupancy of B sites with spins in the state 0, Pb(0). 

An example for such anomalies is the overshooting of \t71b\(T) at low temperatures 
in the vicinity of the 01 ground state, illustrated in figures 2 and 3. As discussed above, 
the overshooting is related to decreasing the number of spins in state 0, by flipping such 
spins to, say, —1, costing minimal energy. Indeed, the anomaly in |ms| is monitored to 
be accompanied by a pronounced decrease in p_e(0). 

Anomalies may also occur in the specific heat, being caused either by a rapid 
increase or decrease in the number of S=l spins in state 0. For instance, such an 
anomaly in C has been observed before [10] in the AF phase at J2/J1 = 0, close to the 
Oil structure. There, monitoring C(T), a three-peak-structure of the specific heat has 
been found, with a non-critical maximum at low temperatures, due to easy flips of spins 
in the state 0, followed, at higher temperature, by a strongly size dependent critical 
peak, and, at even larger temperature, another non-critical maximum, due to flipping 
of single spins for sublattices with rather large clusters of '+' or '— ' spins. 

We observe a similar thermal behaviour of the specific heat in the 01 phase, when 
turning on the nnn coupling, J 2 < 0. An example is depicted in figure 9, at J2/J1 = —0.2 
and D / J\= 3.4. However, here the peak at lowest temperature, Ti, is due to a rather 
drastic decrease of B spins in state 0, which has been discussed above. Actually, by 
enhancing the strength of the single-site anisotropy term, D/Ji = 3.6, the position of 
that peak, at 7], shifts first to somewhat higher temperatures. It still corresponds to a 
lowering, with temperature, of the average number of 0's for the sublattice B, p B (0)(T). 
When moving even closer to the Oil structure, D / J\ = 3.8 and 3.9, 7] tends to be shifted 
towards lower temperatures with increasing D. The peak position seems to follow the 
same dependence as for vanishing J 2 pH] , namely k bTi / J\ oc 4 — D / J\. As may be seen 
from monitoring Pb(0), the peak is then, as in the limit J2 = 0, due to an increase in 
the number of B spins in state at low temperatures. 

5. Summary 

We have studied a mixed spin Ising model with antiferromagnetic couplings, J\, between 
spins S= 1/2 and S= 1 on neighbouring sites of a square lattice, augmented by couplings, 
J2, between spins S= 1 on next-nearest neighbouring sites of the lattice. An additional 
quadratic single-site anisotropy term, D, acts upon the S=l spins. Based mainly on 
ground state considerations and on extensive standard Monte Carlo simulations, we have 
determined the ground state phase diagram in the (D/Ji, J2/J1) plane and identified 
compensation points, types of phase transitions corresponding to different ground state 
structures as well as anomalies for various physical quantities. 

In particular, compensation points are found to exist, in contrast to previous belief. 
They exist for antiferromagnetic nnn couplings, J2 > —1.0, in the 01 phase, springing, at 
zero temperature, from the line J2/J1 = —2 + 79/(2 Ji). Across that line, different types 
of one-spin flips on the S=l sublattice cost lowest energy, accompanied by a change in 
the sign of the slope of the magnetization of that sublattice at low temperatures. At fixed 
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value of J2/J1, the range, in D/Ji, of compensation points is quite narrow, becoming 
extremely small for decreasing nnn couplings, as we have shown by decreasing J2/J1 
from -0.2 to -0.8. 

The transition to the paramagnetic phase appears to be always in the universality 
class of the two-dimensional Ising model, with the critical maxima of the specific 
heat growing with lattice size, L, in a logarithmic fashion, and with the maxima 
of the appropriate (staggered) sublattice susceptibilities growing with lattice size 
proportionally to L 7 / 4 . These results hold in case of AF, 01, and B-AF ground states of 
the model. In case of the Oil ground state, there is no phase transition. 

The model is observed to exhibit interesting anomalies by showing non-monotonic 
temperature dependences for various quantities, including sublattice magnetization and 
specific heat. The anomalies are typically caused by energetically favoured flips of S= 
1 spins between the state and a non-zero state. 

In conclusion, our study provides clear evidence for non-expected compensation 
points in a rather simple mixed Ising model with antiferromagnetic next-nearest 
neighbour couplings between S=l spins on a square lattice. Extensions to three 
dimensional lattices are well beyond the present study, and may be investigated in 
the future. 
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